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We propose the viscous Camassa-Holm equations as a closure approximation for the Reynolds- 
averaged equations of the incompressible Navier-Stokes fluid. This approximation is tested on tur- 
bulent channel flows with steady mean. Analytical solutions for the mean velocity and the Reynolds 
shear stress across the entire channel are obtained, showing good agreement with experimental mea- 
surements and direct numerical simulations. As Reynolds number varies, these analytical mean 
velocity profiles form a family of curves whose envelopes are shown to have either power law, or 
logarithmic behavior, depending on the choice of drag law. 



Laminar Poiscuillc flow occurs when a fluid in a 
straight channel, or duct, is driven by a constant up- 
stream pressure gradient, yielding a parabolic streamwise 
velocity profile which is symmetric about the midplane 
of the channel. In turbulent states, the mean streamwise 
velocity profile remains symmetric, but is flattened in the 
center because of the increase of the velocity fluctuation. 
Although a lot of research has been carried out for tur- 
bulent channel flow, e.g., accurate measurement of 
the mean velocity and the Reynolds stress profiles, in 
particular for flows at high Reynolds numbers, is still an 
experimental challenge and the fundamental understand- 
ing of how these profiles change as functions of Reynolds 
number is still missing. 

In turbulent channel flows it is customary to define a 
characteristic velocity it, and Reynolds number Rq by 
u* = \/\ T a\l P and Rq — du^/v, where To is the bound- 
ary shear stress, we take the density p to be unity, d is 
the channel half-width and v is the molecular viscosity 
of the fluid. Based on experimental observation and nu- 
merical simulation, a piecewise expression of the mean 
velocity across channel has been proposed, for which the 
nondimensional mean streamwise velocity, <f> = U/u„, is 
assumed to depend on 77 = u^z/v and have three types of 
behavior depending on the distance away from the wall 
boundary, z: a viscous sublayer, in which <j) ~ 77; the von 
Karman-Prandtl logarithmic "law of the wall," in which 
(f>(rj) = K _1 lnr7 + A, where K ~ 0.41 and A ~ 5.5; and a 
power law region, in which <f> ~ if, 0<p<l§. 

In this paper, we propose the viscous Camassa-Holm 
equations (VCHE) in (||) as a closure approximation for 
the Reynolds equations. The analytical solutions of the 
steady VCHE depend on the three parameters: the flux 
Reynolds number R — dujv (where u is the streamwise 
velocity, averaged across the channel); the drag law for 
wall friction D = 2tq/u 2 — 2R\jR 2 as a function of R; 
and a shape factor c G (0, 1) which specifies the flat- 
tening of the velocity profile and is given by the ratio 



Umax/u = (3 — c)/2 of the maximum streamwise velocity 
Umax (occurring at the center of the channel) to its av- 
erage u. This family of analytical solutions of the steady 
VCHE provides profiles of the mean velocity and the 
Reynolds shear stress depending on the three parameters 
R, D and c. (Giving these three parameters is equivalent 
to giving the Reynolds numbers R, Rq and R m ax, cor- 
responding to the channel flux, boundary stress and ve- 
locity at the midplane, respectively.) The VCHE profiles 
agree well with data obtained from turbulent channel flow 
measurements and simulations across the entire channel. 
As the Reynolds number is varied in these solutions, the 
velocity profiles <p(rj, R) form a family of curves with up- 
per and lower envelopes. The Blasius i? -1 / 4 law [|| for 
the drag coefficient of the wall boundary leads to rj 1 / 7 
power law behavior for these envelopes, while the von 
Karman drag law leads to logarithmic envelopes. 

We begin our theoretical treatment by recalling the 
Reynolds-averaged Navier-Stokes equations 0j8| 

B{u) 



dt 



(u) • V(u) = div (T), div (u) = 0, (1) 



where (u) denotes the ensemble average of the veloc- 
ity and (T) = -(p)I - (u <g> u) + z/(V(u) + V(u) T ). 
For turbulent channel flow, the mean velocity is of the 
form (u) = (U(z), 0, 0) T , with (p) — P(x,y,z) and the 
Reynolds equations ([!]) reduce to div (T) = 0, or equiva- 
lently, 

— vU" + d z (wu) — —d x P, 
d z (wv) = -d y P, d z (w 2 ) = -d z P, (2) 

where (u, v, w) T is the fluctuation velocity in the infi- 
nite channel {(x,y,z) G K 3 , —d<z< d}. The (1,3) 
component of the averaged stress tensor (T) is given by 
(T13) = —vU{z) + (wu). On the boundary, the velocity 
components all vanish and one has the stress condition 



T T = (T13) 



z = ±d 



vXJ (z)\ 



=±d ■ 



(3) 



1 



upon using (wu) — at z — ±d. Hence, the Reynolds 
equations imply (wv)(z) = and P — Pq — T$x/d — 
(w 2 )(z), with integration constant Pq. 

The viscous Camassa-Holm equations (VCHE) are 



dv 
~dt 



+ VjVu 3 + V (p 



1 



|Vu| : 



(4) 



with V • v = = V • u, where v = (1 — a 2 A)u, 
d/dt = d/dt + u • V is the material derivative. Also, a is 
a constant lengthscale, |Vu| 2 = Tr(Vu ■ Vu T ) and A is 
the Laplacian. Equation (Q) with v = is derived in || 
by decomposing Lagrangian parcel trajectories into mean 
and fluctuating parts, then applying asymptotic expan- 
sions, Lagrangian means, and assuming isotropy of fluc- 
tuations in Hamilton's principle for an ideal incompress- 
ible fluid. That derivation generalizes a one-dimensional 
integrable dispersive shallow water model studied in Jict l 
to n-dimensions and provides the interpretation of a as 
the typical amplitude of the rapid fluctuations over whose 
phase the Lagrangian mean is taken in Hamilton's prin- 
ciple. Moreover, that derivation makes it clear the solu- 
tions of VCHE are mean quantities. 

Before comparing VCHE with Reynolds averaged 
equations, we rewrite (g) in the equivalent 'constitutive' 
form 



du 

~dt 



div T, T = -pi + 2v(l - a 2 A)B + 2a 2 D, (5) 



with Vu = 0, D = (l/2)(Vu + Vu T ), ft = (l/2)(Vu- 
Vu T ), and co-rotational (Jaumann) derivative given by 
D = dD/dt + Dft — flD. In this form, one recognizes 
the constitutive relation for VCHE as a variant of the 
rate-dependent incompressible homogeneous fluid of sec- 
ond grade [jllj, [^2| , whose viscous dissipation, however, 
is modified by the Hclmholtz operator (1 — a 2 A). There 
is a tradition at least since Rivlin (l|] of modeling tur- 
bulence by using continuum mechanics principles such as 
objectivity and material frame indifference. For example, 
this sort of approach is taken in deriving Reynolds stress 
algebraic equation models fuj . Rate-dependent closure 
models of mean turbulence have also been obtained by 
the two-scale DIA approach Jl5| and by the renormal- 
ization group method Since VCHE describe mean 

quantities, we propose to use (||), or equivalently (g), as a 
turbulence closure model and test this ansatz by applying 
it to turbulent channel flow. The corresponding closure 
for turbulent pipe flows will be treated elsewhere jnj . 

We denote the velocity u in (|J) by U and seek its 
steady state solutions in the form U = (U(z),0, 0) T sub- 
ject to the boundary condition U(±d) = and the sym- 
metry condition U(z) — U{—z). In this particular case, 
the steady VCHE reduces to, 



-vU" + va 2 U"" = -d x p, = 
= -d z (p-a 2 (U'f 



(6) 



That is, dv/dt vanishes in (Q) and the remaining a 2 terms 
modify the pressure and dissipation. Comparing (^) and 
(p), we identify counterparts as, 



U — U, d z (wu) — va 2 U + Pq , d z (wv) = . 
V(P+(u- 2 » = V{p-p x-a 2 {U') 2 ), 

for a constant po- This identification gives 

(wv)(z) = 0, —(wu)(z) = —poz — ua 2 U"' (z) , 



(7) 



(8) 



and leaves (w 2 ) undetermined up to an arbitrary function 
of z. A closure relation for —(wu) involving U"'(z) also 
appears in Yoshizawa jL5|, cf. also equation (8) of Q. 

The solution of the steady VCHE (g) in a channel sub- 
ject to these boundary and symmetry conditions |18| is 



U{z) = a(l 



cosh(z/a) 
cosh(d/a) 



6 1 



z 
rf 2 



(9) 



with constants a, b. Any time dependent solution 
of (g) such that u(z,t) = (U(z, t), 0, 0) T with U(z,t) = 
U(—z,t) and U(±d, t) = converges exponentially in 
time to the solution (g) with the same mean flow and 
boundary shear stresses jnj . 

We consider channel flows for R >> 1. Our basic as- 
sumption is that £ = d/a — ->ooasi?^oo. The solution 
@ must satisfy the stress condition (g), which imposes 



TTl, M aV 1 J- 

t = vU {z) \ z =-d = — tanh£ + — — . 

a d 



Substituting the definitions 

rd 



i y"^[/(z)dz = a(l-r 1 tanhe) + 



(10) 



(11) 



and R 



ud 
v 



Ra — 



r n a 



2bv 

dr a ' 



with < 6 < 1, into relation (|l^) re-expresses it as 



Rq = (R — 



ORf 



^ tanh £ 



1 - ^tanh^ 



QR n 



(12) 



Upon setting ^ = SRq and taking ^ >> 1, this simplifies 
to order 0(l/£) into the basic relation 



1 - . 



R @Rq 

Rq 3 



(13) 



In terms of the parameters 8 and £, equations (@) and (j|) 
imply the Reynolds shear stress, 



- (wu)(z) = T (I - I 



sinh(z/a) 



sinh(^) 



(14) 



Thus, the solution (ph implies —(wu) > for —d < z < 0, 
as seen empirically [H . In the lower half of the channel, 



2 



the symmetric solution ([)[) may be expressed in wall units 
using the notation (f>(rj) — U(z)/u*, rj — (z + d)/£*, with 
£* = vju* = d/Ro, to order 0(£e~£) as 



077 1- 



2Rc 



(15) 



for < rj < Rq. In this notation, we have a = d/£ = £*/5 
for the lengthscale in (^). The velocity profile U{z) in (||) 
has its maximum at the center of the channel (rj = Rq). 
At this point <f> max — 4>{Ro) is given to leading order by 



1 - i 



QRp 
2 ' 



(16) 



Hence, from (gj), (|§) and ^/2/D = R/R we have 
6Rq = 6( 



2 ^ 1 ~ d W 2 <M 



(17) 



Since < 6 < 1, relations (|T^) imply the inequalities 
3/2 > (f>max\/ D/2 > 1, and we may write 



/2 1 
— , and — 



(18) 



by introducing the quantity c £ (0,1) defined in terms 
of the velocity profile flatness or midplane velocity ratio 

Umax I H &S 



3 - c 



■u,a.r.\l— With 0<C<1. (19) 



Comparison with data will show that c is nearly indepen- 
dent of R and bounded away from its endpoint values. 
The first equation in (|lj) and the basic relation (|ll|) then 
imply 



O(i). 



(20) 



3(1 -c) 

Substituting this into the second equation in (111) gives 



i?o = c6R(i 



3 -^)=c5R + 0( 1 -). (21) 



Thus, to leading order, (5 = c 1 Rq/R = c 1 y/D/2 and 
the velocity profile ([l?]) is given by 



R 
Rq 



!(l _ e -* )+3(1 -^ (l -i) 



(22) 



Thus, for £ = cZ/a >> 1, the drag y/D/2 = R /R and 
the constant c determine the steady velocity profile of 
VCHE ip(rj) at each R. The lengthscale a is given to 
leading order by a/d = 1/(5R ) = 2c/ \DR). For the 
Blasius drag law, D — Ai? -1 / 4 , A = const, this implies 



H=f- 3/4 =<!) 4, V" 



(23) 



This scaling in R agrees with the scaling law for the Kol- 
mogorov fluctuation dissipation length, £d- Thus, in this 
case, a is proportional to £ d H]. For the Blasius drag 
law, we also have from (n8h 



_ 3(1 - c) 3/4 1-0 _ c 
A/2 ' 5 ^\/2 



R L 



(24) 
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FIG. 1. The mean- velocity profile <j> as a function of log(r7). 
Experimental data |Q| , DNS data |J . (f> v and 0t represent the 
viscous relation and von Karman logarithm law respectively. 
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FIG. 2. The Reynolds-stress — (iou)/u». Experimental 
data [@, DNS data §. 

Figures 1-3 compare our theoretical solution ( p2] ) with 
experimental data from [0, and direct numerical simula- 
tion (DNS) data from [pL We assume the Blasius drag 
law, D = Ai? -1 / 4 , with A = 0.06. We vary c slightly with 
R in order to best fit the data (c S [.728, .77]). Note, 



3 



the curves can be brought into even better agreement by 
allowing A to also vary slightly with R. 
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FIG. 3. The turbulent kinetic energy production profile. 
Experimental data is from Wei & Willmarth W. 
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FIG. 4. The upper and lower envelopes of the velocity pro- 
file for the Blasius drag law. 

After the standard Blasius drag law is chosen and c is 
determined from the midplane velocity data, no free pa- 
rameters remain in the model. Figure 2 compares the 
theoretical prediction with the measured data for the 
Reynolds-stress ~(wu)/ul. Figure 3 compares the pre- 
dicted and measured turbulent kinetic energy produc- 
tion profiles. Except for the region nearest the wall, 
log 10 (?/) < 1 in Figure 2 and Figure 3, the theoretical 
model shows excellent agreement with the data. Using 
the fundamental relation (^J) yields a family of velocity 
profiles 4>(r/;R), plotted in Figure 4 at various values of 
R for the Blasius drag law. Remarkably, the upper and 
lower envelopes of this family are found analytically to be 
77 1 / 7 power laws. If instead of the Blasius drag law, one 



uses the von Karman drag law, \[2/D = Ai log(i?o) + A 2 , 
for constants Ai and A2, then (|2^) yields a family of ve- 
locity profiles 0(?7; R) whose upper and lower envelopes 
are nearly linear in log 10 7/ with different slopes. 

We are grateful to R. Kraichnan for constructive com- 
ments and to D. Cioranescu for pointing out the relation 
between the Camassa-Holm equations and the second- 
grade fluids. 
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